Effect of pinning on the vortex-lattice melting line in type-II superconductors 
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The vortex-lattice melting line in three-dimensional type-II superconductors with pinning is de- 
rived by equating the free energies of the vortex system in the solid and liquid phases. We account 
for the elastic and pinning energies and the entropy change that originates from the disappearance 
of the phonon shear modes in the liquid. The pinning is assumed to be caused by point defects 
and to be not too strong so that the melting line lies inside the so-called bundle-pinning region. 
We show that the derived equation for the melting line is equivalent to some Lindemann criterion, 
which however differs from that used previously. Estimating the effect of pinning on the entropy 
jump at melting, we find the upper critical point of the melting line from the condition that this 
jump vanishes. We also consider the H-T phase diagrams of type-II superconductors for different 
strengths and types of pinning and analyze the two recently discussed scenarios how the melting 
line and the order-disorder line merge. 

PACS numbers: 74.25.Qt, 74.72.Bk 
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I. INTRODUCTION 

In three-dimensional high-T c superconductors with 
pinning, two phase transition lines are known to exist 
in the magnetic field H - temperature T planei 2 ^*^ 
The line H m (T) where a quasiordered Bragg glass&£ ther- 
mally melts into a flux-line liquid, and the order-disorder 
transition line Hdi S (T) separating the Bragg glass from 
an amorphous vortex state. The melting is caused by 
thermal vibrations of the lattice, while the order-disorder 
transition is induced by quenched disorder in the vortex 
system. These two lines merge at some point in the H-T 
plane. Although both transitions are accompanied by a 
proliferation of dislocations in the vortex lattice, it was 
argued^ that the dislocation density p is essentially dif- 
ferent in these cases: p ~ a~ 2 for melting, and p ~ R~ 2 
for the order-disorder transition. Here a is the spac- 
ing between flux lines, and R a is the so-called positional 
correlation length 9 within which the relative vortex dis- 
placements caused by the quenched disorder are of the 
order of a. In fact, an intersection of these two differ- 
ent phase transition lines occurs in this scenario, and the 
order-disorder line terminates at the intersection point 
while the melting line continues for some distance to 
higher H, see Fig. ^ Within this physical picture, the 
existence of the so-called slush phased can be naturally 
explained. Recent experimentsAi*i2iiii^ii^ for YBaCuO 
seem to support this scenario. On the other hand, exper- 
imental dataiSiiLiS for BSCCO strongly argue in favor 
of a diffe rent scenario which was implied, e.g., in the 
Refs. I2l5i In this second scenario, the dislocation densi- 
ties for both lines coincide at the point where these lines 
merge, and in fact, one deals with only one phase transi- 
tion line that describes both the order-disorder transition 
at low temperature and the melting near T c , Fig. ^ 

Phase diagrams of superconductors with pinning re- 
flect the competition of three characteristic energies^: the 



elastic energy, the pinning energy, and the energy of ther- 
mal fluctuations. At melting, the cost in the elastic en- 
ergy due to the proliferation of dislocations is mainly bal- 
anced by the entropy gain associated with thermal fluc- 
tuations, while the role of the pinning energy, as shown 
below, is determined by the parameter a/R c . Here R c is 
the transverse collective pinning length^. On the other 
hand, at the order-disorder transition, the balance of pin- 
ning energy and elastic energy is most important, while 
the relative contribution of the entropy gain is negligi- 
ble at low temperatures and, according to the scenario of 
Ref.|8|, is determined by the ratio a/R a near the intersec- 
tion point. Thus, if the intersection of the melting and 
the order-disorder lines occurred sufficiently deep in the 
bundle pinning region (so that R c 3> a at this point), the 
scenario of Ref. |3 would lead to the conclusion that up 
to the intersection point, one can find the melting line 
by neglecting the pinning, and the order-disorder line 
by neglecting the entropy gain. Just this approximation 
was used in our paper— for analyzing the phase diagrams 
of superconductors. However, our recent results^ point 
out that flux-line pinning can affect the melting line near 
the intersection point since at this point, the ratio R c /a 
has decreased to several units for any magnitude of the 
quenched disorder in the vortex lattice (even when the 
disorder is weak). As to the second scenario, the three 
energies are all the same order of magnitude in the tem- 
perature region where the order-disorder transition grad- 
ually transforms into melting. Thus, whatever scenario 
occurs in reality, a detailed investigation of the effect of 
pinning on the melting line is important to clarify the 
most intriguing part of the phase diagram. 

The effect of pinning by point defects on the 
melting line was observed both in BSCCO 20 and in 
YBaCuQ2i*22*2£ crystals. It was discovered that an in- 
crease of the quenched disorder in the vortex lattice leads 
to a noticeable shift of the intersection point to lower 



2 



magnetic fields and simultaneously pushes the melting 
line in the T-H plane slightly downwards, i.e., at a fixed 
temperature the appropriate magnetic field of the melt- 
ing, H m (T), decreases. It is important that the shift 
of the intersection point is essentially more pronounced 
than the decrease of H m {t) itself. 

Some theoretical results on this subject were obtained 
in Refs. I2ll5ll24l Larkin and Vinokur— started from the 
assumption that for the vortex lattice with quenched dis- 
order to melt, the temperature must match a character- 
istic barrier composed of the elastic energy, E e \, and the 
pinning energy, E p - ln . So they estimated the effect of 
pinning on the melting line by considering the following 
balance of these three energies: 

T = E pin + E c \, (1) 

where E pm and E e \ were calculated in the so-called cage 
modelji But it follows from this equation that the melt- 
ing line has to shift upwards when the quenched disorder 
increases. Another approach was used in Ref. y. To de- 
scribe the unified phase transition line, Giamarchi and Le 
Doussal, 2 who implied the second scenario, put forward 
the following generalization of the Lindemann criterion: 

2 2 2 m\ 

"toted = C -L° J ( 2 ) 

where cl is the Lindemann constant (cj, ~ 0.1 — 0.2), 
and utotai ~ [211^ + u 2 (a,0)] 1 / 2 is the rms displacement 
of neighboring flux lines caused both by the thermal fluc- 
tuations and by the quenched disorder in the lattice. 
Here ut is the magnitude of the thermal fluctuations, 
while u(a, 0) describes the mean relative displacement of 
neighboring flux lines caused by the disorder. This cri- 
terion leads to the usual condition for the order-disorder 
transitioni2*4 at low temperatures when ut is negligible, 
and it goes over to the well known Lindemann criterion 
for pure melting when the disorder disappears. Equa- 
tion J2| results in a qualitatively correct dependence of 
the melting line on pinning by point defects. However, 
this reasonable criterion is only an interpolation formula 
between the two limiting cases and has no serious jus- 
tification. The same is true for the criterion of Ref. 5. 
Radzyner et al£ described the unified phase transition 
line using the criterion, 

u 2 T + u 2 {a,0) = c 2 L a 2 , (3) 

which practically coincides with Eq. (J5J. The criterion 
© is equivalent to the following balance of energies: 

T + E pin = E ch (4) 

which evidently differs from Eq. ft}. In calculations of 
the phase transition line, the energies -E p i n and E c \ were 
estimated in Ref. |B| in the framework of the cage model. 
But it remained unclear why the pinning energy E pm 
now enters into Eq. (0J with the opposite sign as com- 
pared to the energy balance (J2J suggested by Larkin and 
Vinokur iw 
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FIG. 1: Schematic T-H phase diagrams for the first (solid 
lines) and for the second (dashed line) scenarios. In the first 
scenario the melting line terminates at the so-called upper 
critical point (T up , H up ) which in general does not coincide 
with the intersection point (T, Hi). In this case a slush phase 
(i.e., a vortex liquid with smaller density of dislocations) can 
be observed. In the second scenario the order-disorder and 
the melting lines are manifestations of a unified phase transi- 
tion line. Note that for both scenarios the vortex liquid and 
the amorphous vortex state are, in fact, one and the same 
phase, which has different viscosity at low and high tempera- 
tures. (We do not discuss the "vortex glass transition" which 
may be not a true phase transition.) 



In the present paper, in order to find and to justify 
a criterion for the melting of the flux-line lattice with 
quenched disorder, we start with the analysis of melt- 
ing in the ideal lattice and show that three different ap- 
proaches lead to the same dependence H m (T). These 
approaches are: the Lindemann criterion, the energy bal- 
ance, and the rigorous approach based on the Ginzburg- 
Landau (or on the London) Hamiltonian. In the case of 
the lattice with quenched disorder, we show that the pin- 
ning energy in the flux-line liquid is larger than the pin- 
ning energy in the Bragg glass. For this reason, the dif- 
ference of these pinning energies, E p i n , has the opposite 
sign as compared to Eq. ft}, and we arrive at an equation 
similar to Eq. but with an expression for E p i n that 
differs from the estimates previously published. Besides 
this, we estimate the effect of pinning on the entropy gain 
at melting and find the upper critical point of the melting 
line from the condition that this gain vanishes. We also 
show that the result for H m (T) based on this energy bal- 
ance agrees with the result which can be derived from the 
Ginzburg-Landau Hamiltonian. Then, using the H m (T) 
obtained in the framework of the second and the third ap- 
proaches, we find how the Lindemann criterion should be 
modified to give the same melting line. Interestingly, the 
presented energy balance clarifies the difference between 
the first and the second scenarios. Finally, we present 
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the T - H phase diagrams of superconductors with two 
types of fiux-line pinning by point defects and compare 
these diagrams for the two scenarios. 

In this paper we consider only magnetic fields exceed- 
ing considerably the lower critical field H c i and thus do 
not distinguish between the magnetic field H and the 
magnetic induction B. Besides this, we deal only with 
uniaxial anisotropic three-dimensional superconductors, 
neglecting completely the decoupling of the supercon- 
ducting layers. The anisotropy is characterized by the 
parameter e = A a b/A c < 1 where X a b and A c are the 
London penetration depth in the plane ab perpendicular 
to the anisotropy axis and along this axis, respectively. 
The magnetic field is assumed to be directed along the 
anisotropy axis. As to the quenched disorder in the flux- 
line lattice, we assume that it is caused by point defects 
and is not too strong such that the melting line lies en- 
tirely in the bundle pinning region. 



II. MELTING OF THE IDEAL VORTEX 
LATTICE 

We begin with the analysis of melting in the ideal 
pinning-free vortex lattice and compare the results of var- 
ious approaches. 



A. Lindemann criterion 

According to the well known Lindemann criterion, the 
flux-line lattice melts when the magnitude of the thermal 
displacements of the lattice relative to its equilibrium 
position, ut, reaches a certain fraction of the spacing 
between the flux lines, o: 

2 2 2 rr\ 

u T = c L a , (5) 

where a = ($0/ H) 1 / 2 , $0 is the flux quantum, and cl 
is the Lindemann constant. The magnitude ut depends 
on the elastic moduli of the latticed and was calculated 
in many papers; see, e.g., Refs. 0, |2f|26 27 .28. It can be 
represented in the formi* 



Gi 



l-t 2 



1/2 



h- 1/2 f(h), 



(6) 



where £(t) is the coherence length in the ab plane, h = 
H/H c2 (t), t = T/T c , H c2 (t) = $ /2< 2 is the upper 
critical field, Gi is the Ginzburg number, 



Gi 



1 / T c 



2 V#c 2 ^o 3 



which characterizes the width of the fluctuation region in 
zero magnetic field, £0 and H c are the coherence length 
and thermodynamic critical magnetic field of the super- 
conductor in the Ginzburg - Landau theory extrapolated 
to T = 0. For dcfmitcness, we implied in Eq. © and 



below that £ 2 (t) = £ 2 (0)/(l - t 2 ). [Hence £ = £(0)/a/2]. 
The complete expression for the function f(h) was given 
in Ref . I28L but for our further analysis it is sufficient to 
use a simplified form of this function* in which the contri- 
bution containing the compression modulus of the vortex 
lattice, en, is neglected: 



f(h) = 



2/3 A [l + (l + 5) 2 ] 1/2 -l 
1 - h c(l + c) 



(7) 



with 5 = 0.5[/3 A (l - h)} 1 / 2 , and A = 1.16. Note that 
this formula can be rewritten in the form: 



f(h) = 



{1-hf/ 2 '- 



(8) 



where the function f\{h) defined by this equality de- 
creases monotonically with increasing h, and its varia- 
tion in the interval < h < 1 is not large: /i(0) « 2.34, 
/i(l) ~ 1.78. Thus, to a first approximation, this func- 
tion can be considered as a constant, fi(h) w /1 ~ 2. 

Combining formulas J^J-Q, we arrive at the equation 
for the normalized melting field h m (t) — H m (t) / H c2 (t) : 



t 



Gi 
l-t 2 



1/2 



This equation agrees with those obtained 
earlier^^SiSL 2 ^ and differs from them only in the 
form of the function fi(h) since different authors used 
slightly different approximations for the elastic moduli 
or took into account the contribution associated with 
the compression modulus of the vortex lattice. 

When the normalized melting field h m (t) is small, 
h m <C 1, it follows from Eq. 10 that 



H m (t) 

H c2 {0) : V/i(0)f 



27TC? \ 2 (1 - t 2 ) 2 



Gi 



(10) 



This is a well-known result 9 i 25 i 26 i 27 ' 28 that holds for tem- 
peratures near T c (but outside the fluctuation region in 
zero magnetic field), namely for 



Gi < 1 - t 2 < (/ 1 (0)/27rc|) 2 Gi. 



(11) 



Note that for inequalities l|llfl to be fulfilled in a suffi- 
ciently wide temperature interval, the Ginzburg number 
should not be too small. In the opposite limiting case, 



l-^» (f 1 (l)/2nc 2 L ) 2 Gi 



(12) 



the field h m is close to unity, and one obtains from 
Eq. ©: 



1 



or equivalently, 



{2nc 2 L 



2/3 



.2/3 



Gi 
l-t 2 



1/3 



(13) 



H c2 (t)-H m (t) //i(l)\ 2/3 , 



H c2 (0) 



V27TC 2 



t 2/3 G .l/3 (l _ t2) 2/3 (u) 
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Note that H c2 (0)Gi 1/3 (l-t 2 ) 2 / 3 is the width (along the H 
axis) of the fluctuation region in not too small magnetic 
fieldsj 29 ! 30 ! 31 H > GiH c2 (0), and for inequality O to 
hold, the Ginzburg number should not be too large. 



B. Energy balance 

At melting, proliferation of dislocations occurs in the 
vortex lattice. These dislocations create a network in the 
lattice, and we consider a mean unit cell of this network 
composed of edge and screw dislocations. The energies 
of these dislocations are of the order of 



-Eodgo ~ C66 CI 2 1 1| , 
-Escrow ~ (C44C66) 1 ^ 2 2 ^_L, 



(15) 



where Cqq and C44 are the shear and tilt moduli of the 
flux-line lattice, while i|| and l± are the dimensions of 
the cell in the longitudinal and transverse directions to 
H, respectively. Since at the melting these dimensions 
are of the order of a (see below), we have omitted the 
logarithmic factors ]n(lu/a), hi(l±/a) in the above for- 
mulas for the dislocation energies. The shear and tilt 
moduli may be expressed a,sfi^ 



C66 p 

c 44 



e 2 e 



(1-h), 



where h = H/H c2 (t), £0 = (<E>o/47rA a fc) 2 , and we have 
omitted the logarithmic factor of type ln(a/£) in the tilt 
modulus (this nonlocal modulus should be estimated at 
wave vectors k of the order of 17 , /J 1 ). The factors 
containing (1 — h) take into account the softening of 
the vortex lattice near the H C 2(t) linei 3 ^ Minimization 
of the elastic energy of the cell at its fixed volume leads 
to -B^dgc ~ -Escrow and hence gives the relation between 
Z|| and l±: 



C44 

C66 



1/2 



(17) 



Thus, up to a numerical factor, the elastic energy of the 
dislocation cell, E c \, equals E SCICW . 

At melting, the cost in the elastic energy due to the 
proliferation of the dislocation network is balanced by 
the entropy gain in the free energy of the flux-line lattice. 
We now give a simple estimate of this gain: The vortex- 
lattice degrees of freedom associated with shear undergo 
a change when the melting occurs. This change occurs 
for lattice modes with wave-lengths greater than l±. 
There are modes of this type in the unit volume of 

the lattice, and each of them contributes about T to the 
entropy gain. Thus, up to a numerical factor, the gain 
per cell of the dislocation network is T -ZiT^T 2 ~ T. 



We now can write down the change of the free energy 
per cell at melting: 



AF cx C{c 44 c 66 ) 1/2 a 2 l ± - T, 



(18) 



where some constant C is the ratio of the unknown nu- 
merical factors in E e \ and in the entropy gain. Minimiza- 
tion of Eq. l|18|l with respect to the parameter l± leads to 
the conclusion that this parameter should have the min- 
imum possible value. It is clear that this value is of the 
order of a in the lattice. Then, taking into account that 
AF = at the melting, we arrive at 



C(c44C 6 6) 1/2 a 3 - T = 0, 



(19) 



where this constant C may be slightly renormalized as 
compared to the C of Eq. (|18fl . Inserting Eqs. (|16fl for 
the elastic moduli and assuming X a b(t) = const, one 
finds the equation for the melting line h m (t): 



1/2 



1/2 



Gi 



1-t 2 ) (I-M 3 / 2 



(20) 



As was noted earlier (see, e.g., Refs.0,0), formula (jTU|) 
for the melting line near T c can be obtained from the en- 
ergy balance. Here we have taken into proper account 
the softening of the elastic moduli near H c2 (T), and now 
equation (|20|l shows that not only formula (|10fl but also 
expression (|13fl can be derived by this method. Moreover, 
if one uses the approximation in which the function fi(h) 
(16) i s a constant^ 3 - fx(h) ~ fi, and put C = 8^c 2 L /fi, 
Eqs. © and (|20() completely coincide in the whole tem- 
perature interval. 



C. Some exact results 

Within the mean-field theory, when one neglects fluc- 
tuations of the superconducting order parameter, the 
melting line H m (T) coincides with the H C 2(t) line. It is 
the fluctuations that shift H m (t) downwards in the H-T 
plane. As was mentioned above, see Eq. (|14[l. the Lin- 
demann criterion shows that at sufficiently strong mag- 
netic fields, the distance between the melting line and the 
mean- field H c2 (t) line is comparable with the width of the 
fluctuation region. But then the question arises about the 
applicability of this criterion (and of the energy balance) 
for determining H m in this region of the magnetic fields 
since expressions Ijl6|l for the elastic moduli were derived 
in the framework of the mean-field theory without ac- 
counting for the fluctuations. In this context, it should 
be noted that one cannot confine oneself to taking into 
account only the first fluctuation correction to the elas- 
tic moduli (in the fluctuation amplitude) since inside the 
fluctuation region the amplitude is large, and corrections 
of all orders are essential. In particular, the renormalized 
elastic moduli will vanish on a line which differs from the 
mean- field H c2 (t). However, simple considerations 3 ^ 3 ^ 
show that at 1 — h m <C 1, strong fluctuations can only 
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renormalize the numerical factor in Eq. (|I3[1 . but the de- 
pendences of H m on t and on Gi remain unchanged. We 
now briefly outline these considerations. 

In dimensionless units the Ginzburg-Landau Hamilto- 
nian depends on the three parameters: t, H/H C 2(0) and 
the Ginzburg number Gi. As well known, the quadratic 
part (in the order parameter) of this Hamiltonian looks 
like the Hamiltonian of a particle with double electron 
charge in a magnetic field, and the energy spectrum of 
this particle is the so-called Landau levels. It is essential 
that in fields H 3> Gi H C 2(0), the distance between these 
levels exceeds the width of the fluctuation region. Thus, 
if one expands the order parameter into the eigenfunc- 
tions of the particle, only the modes of the order param- 
eter corresponding to the lowest Landau level strongly 
fluctuate near the melting line, and one may retain only 
these modes in the Hamiltonian to calculate the fluctua- 
tion part of the free energy of a superconductor. It turns 
out that the Hamiltonian thus obtained depends on a 
single combination of the parameters. In our notations, 
this combination, Q, can be represented in the form: 



Q = 



(i-fr)(i-< 2 ) 1/3 

£2/3(^1/3^2/3 • 



(21) 



where h = HjH^ii). Hence, in this region of the mag- 
netic fields, the free energies of the vortex liquid, Fu q , 
and the vortex lattice, flat, are also determined only by 
this combination, and for the melting line h m (t) we arrive 
at the equation: 

F liq {Q) = fiat(Q). (22) 
The solution of this equation has the form: 



(i-M(i-Q 1/3 r 



(23) 



where G\ is some constant. Taking into account that 
h rn w 1, we find from equation 123|) : 



1 - h m « dt 2 ^ 



Gi 



1 -t 2 



(24) 



It is seen that up to a numerical factor, this expression 
indeed coincides with Eq. (|13[l . 

The above consideration does not use any perturbation 
theory and therefore is exact, but it does not yield the 
constant C\ . Using some variants of perturbation theory 
(see, e.g., Ref. 135^) . approximate expressions for Fn q and 
-Flat were obtained in Refs. |3(| |33 On this basis, an 
equation for the melting line, which agrees with Eq. J22Jl, 
was derived in these papers together with the appropriate 
constant C\ . Hikami et ali^& estimated C\ ~ 7, while Li 
and Rosenstein^Z. who used a refined expression for F\ at 
found Ci « 9.5. 

A consideration similar to that presented above was 
also applied to the Hamiltonian of the vortex system in 
the London approximation, see Ref. 38. In the region 



of the magnetic fields considered here, H 3> H c %, one 
has X a b 3> a, and the Hamiltonian, as well as the free 
energy of the vortex system, are determined by a single 
combination of the physical parameters: eeoa/T. Thus, 
the melting line H m (T) is found from the equation: 



ee (T)a(H m ) 



= C 2 , 



(25) 



with some constant Ci. In other words, the London 
approximation leads to H m (t) cx [t\ 2 ab {t)]~ 2 ■ Putting 
Aab(i) = A Q b(0)(l — t 2 )^ 1 / 2 and taking into account the 
definition of Gi, we arrive at the formula: 



H m (t) ^ 7T (1-t 2 ) 2 

H c2 (0) {^C 2 t) 2 Gi 



(26) 



which agrees with Eq. I|10|) . 

To summarize, we have shown in this section that up to 
numerical factors, the three different approaches lead to 
the same dependences of the vortex-lattice melting field 
H m on the temperature T and on the Ginzburg number 
Gi, Eqs. IjlOp and (jI4(l . Thus, after obtaining these de- 
pendences, e.g., from the energy balance, for the vortex 
lattice with pinning, one can guess the true form of the 
Lindcmann criterion for this lattice. We shall use this 
procedure in the next section. 



III. MELTING OF THE VORTEX LATTICE 
WITH QUENCHED DISORDER 

In this section we analyze the melting line of the vor- 
tex lattice with pinning assuming that the line is in the 
bundle pinning region (where the transverse collective 
pinning length R c is greater than a). As it follows from 
the figures of Sec. IV B, the melting line, as a rule, does 
entirely lie in this region, and only in the case of suffi- 
ciently strong 5T C pinning can it enter the single vortex 
pinning region. 



A. Energy balance 

We now consider the influence of pinning by point de- 
fects on the energy balance. The adjustment of the vor- 
tex system to the pinning potential decreases the total 
energy of the system, and the amount of this decrease 
is just of the order of the pinning energy. It will be 
shown in this section that at melting, the pinning energy 
in the flux-line liquid, E*^, is noticeably greater than 
the pinning energy in the flux-line lattice (Bragg glass), 
.Epin. Thus, there is a gain in the free energy of the 
vortex system, E p i n , associated with the pinning energy: 



E, 



E^-Elt ~ E^ n . This gain adds to the entropy 



pm -*^pin "^pin ^pin 

gain, and equation l|18|) is modified as follows 



AF cx C(c i4 ,c 6 e) 1/2 a 2 l± - TAS - E } 



pm i 



(27) 



() 



where the factor AS takes into account the effect of pin- 
ning on the entropy gain per dislocation cell. 

In order to estimate E^, and AS", it is neces- 

sary in general to take into account the so-called ther- 
mal depinning 9 '— since the thermal displacement ux is 
sufficiently large at the melting, ut ~ c^a. However, 
to explain the main ideas, we first carry out the esti- 
mates neglecting this dcpinning and then generalize the 
obtained results by taking it into account. 

The pinning energy of the flux-line lattice without dis- 
locations, i.e., of the Bragg glass, in the volume equal to 
the cell of the dislocation network, E^, can be estimated 
using the results of collective pinning theory?2. 



E, 



lat 
pin 



(28) 



where W — f 2 in n^ 2 /a 2 , f pin is the mean elementary 
pinning force exerted by one point defect, n is the con- 
centration of the defects, and the expression (Wl\\l 2 l _) 1 ^ 2 
is the mean pinning force per cell. The displacement 
u = u(lj_, 0) ~ u(0, Zy) is the rms relative shift of two line 
elements in the vortex lattice separated by a distance l± 
transverse to the magnetic field, or by Z|| along the field. 
This shift is caused by the random point defects. The 
magnitude of u can be expressed in terms of the trans- 
verse collective pinning length, i? c , at which the relative 
displacements of points in the lattice are of the order of £. 
In particular, if the small bundle pinning regime occurs, 
i.e., R c < X a b, one h&s& 



e 



1 + HRl/o?) : 



(29) 



where we have used that l± ~ a. Note that according to 
Eq. I|I7|) . the longitudinal dimension l\\ is of the order of 
ea/il-h) 1 / 2 . 

At melting, proliferation of dislocations occurs, and 
in the liquid vortices can adjust themselves to the pin- 
ning potential not only via their elastic deformations as 
in the Bragg glass but also via the plastic vortex-lattice 
displacements generated by dislocations. This additional 
adjustment mechanism increases the pinning energy in 
the liquid. To estimate the pinning energy per dislocation 
cell in the liquid, -E^, it is necessary to take into account 
that the displacements generated by the dislocation net- 
work are essentially larger than the displacements exist- 
ing in the Bragg glass at the same temperature and mag- 
netic field within the scales l±, l^. Indeed, in the lattice 
without dislocations, one has u(l±,0) ~ u(0, In) —u<^ 
[see Eq. (I29[l ]. while in the liquid the dislocations lead to 
displacements of about a > £ at the boundary of the dis- 
location cell. It is also essential that the displacements 
caused by the dislocations within a scale I < l±, l\\ have a 
different dependence on I than the elastic displacements 
in the Bragg glass. Then, we obtain the estimate: 



E 



liq 
pin 



(30) 



Note that although the displacements generated by the 
dislocations are large and exceed £, we multiply the mean 



pinning force per cell only by £ in formula i|3U|) since £ is 
the effective range of the elementary pinning force / p i n - 
It follows from Eqs. pg ]l - PlJ l that E^jE^ n ~ ~ 
[I + ln(i?2/ a 2 )] 1 /2 j i. e . ; ££q noticeably exceeds E p f n in 

the bundle pinning regime, and hence formula (|30|) gives 
the estimate of E pin in Eq. ijTTI) . 

It is convenient to rewrite expression (|30|l in terms of 
L c , the single-vortex collective pinning length, using the 
relation, 9 fp in n e 4 £g/L|?. Then one arrives at the for- 
mula for E-nir,: 



E, 



pin 



ee a[A9o(i)] 3/2 ^ 1/4 (l-M 3/4 , (31) 



where we have inserted the estimates for l±, In, have 
taken into account that the quantity Eq ^ ^ab snou lcl 
lead to an additional factor 1 — h when h = H/H C 2(t) 
tends to unity^ and have used the notation^ 



<m 

L c (t) 



Dg (t) 



(32) 



with e£(0)/L c (0) = D. The function go(t) is given- by 
, 9o (i) = (l-i 2 ) 1/2 (33) 



for SI pinning, and by 

5 o(i) = (l-* 2 )- 1/6 



(34) 



for 5T C pinning. The parameter D is a measure of the 
pinning strength and is estimated as& D i=s (jc/jo) 
where j c is the critical current density in the single vortex 
pinning regime and jo is the depairing current density, 
both taken at T = 0. 

Let us now estimate the entropy term in the energy 
balance (|27fl. In Sec. II B we obtained the expressions 
for the melting line of the ideal vortex lattice, assuming 
that the main contribution to the entropy gain at melting 
is due to disappearance of the shear phonon modes in 
the liquid. The mechanism of this disappearance is the 
following: A shear stress in the liquid generates the so- 
called Peach-Kohler forces^S^i exerted on dislocations; 
the dislocations begin to move, and their shifts relax the 
shear stress in the liquid. In the vortex system with 
pinning a dislocation cannot move if the pinning force per 
dislocation cell, {Wl^l 2 /) 1 / 2 , exceeds the Peach-Kohler 
force acQQU xy l\\ exerted on a dislocation segment of length 
in where u xy is the shear deformation of the vortex system 
in the plane normal to the magnetic field. Thus, for 
relaxation of the shear stress to occur, the deformation 
u xy must exceed the critical value 

= ( W ¥l) 1/2 „ [Dgo{t)] V2 h -m {1 _ h) -m (35 ) 

OC66<|| 

Here we again have inserted the relations used in de- 
riving Eq. I|31|) . In the volume with dimensions L±_ and 
thermal fluctuations generate shear deformations u xy 
that can be estimated from 



c 6e u 



■<!J 



LuL , 



T. 



(36) 
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Hence, for large L± and L«, when 



ULl > [L [{ Ll] 



T 



the deformations are less than the critical value; the re- 
laxation of the shear stress in the liquid does not oc- 
cur, and these modes do not contribute to the entropy 
gain (they do not differ essentially from the appropriate 
modes in the Bragg glass) . Then, applying the consider- 
ations of Sec. II B, we obtain the factor AS" in Eq. J57J): 
AS = 1 - P(t,h), where 



P(t, h) 



hi 2 , 



[Dg Q (t)f{\-t 2 f/ 2 
tGiV 2 h 



(37) 



As is well known thermal fluctuations of the flux- 
line lattice lead to a smoothing of the pinning potential 
and thereby affect the pinning. We take into account 
the effect of this thermal depinning on -E p i n and P(t, h) 
using the recipe of Ref. (see page 1214 in that paper). 
Then an additional factor £/r p appears in Eq. I|3l)|l . while 
the right hand side of Eq. (|35|l is multiplied by (£/r p ) 2 
where r p — (£ 2 + it 2 ") 1 / 2 is the new effective range of the 
pinning force when thermal fluctuations are allowed for. 
This modifies formulas (|31|) and l|37|) as follows: 



E 



pm 



eeoa 



[D 30 (i)] 3/2 ^ 1/4 (l-^) 3/4 



l + t 



Gi 



1/2 



1(h) 



P(t,h) 



l-t 2 J 

D 5o (i)] 3 (l-* 2 ) 1/2 



1/2 



tGi l / 2 h 



l + t 



Gi V /2 f(h)i 2 



l-t 2 



, (38) 



h 1 / 2 



where we have used Eq. © for 

Inserting formulas (|3l3)) in Eq. l|2"7|) . expressing C as 
%^/ttc 2 l / fi with a constant /i, and taking into account 
that AF = at melting, we eventually find the equation 
for the melting line h m (t): 



F T {t)h 



1/2 



(1 - h m )V 2 
A[Dg (t)]y 2 hU 4 



[1-P(t, h m )} 



(39) 



{l-h m f/ 2 + F T (t)h, 



-1/2 



1/2 



2ttc 2 



L ■ 



where 

P(t,h) 



B[Dg (t)] 3 (l-t 2 y/ 2 



tGi 1 ' 2 



h 1 / 2 + F T (t)(l - h)- 3 /< 



(40) 



A, B are some numerical factors, and we have used the 
notation 



F T (t) = fit 



Gi 
~ T 2 



1/2 



(41) 



putting /i = 2 below. Equation generalizes Eq. (|2"U|) 
obtained in the case of the ideal lattice. Note that the 
factor 1 — P = AS in Eq. I|39|) naturally explains the 
existence of the upper critical point on the melting line. 
The temperature t up of this point can be obtained from 
the condition AS" = 0, or 



P(tup> ^m(^up)) 



(42) 



i.e., the entropy jump in the vortex system at melting 
vanishes at this point. 

Equation l|39|) is in agreement with the general con- 
clusion made in Ref. 0. Namely, it was stated^ that the 
phase diagrams of various three dimensional supercon- 
ductors with point defects are determined only by the 
Ginzburg number Gi, the parameter D characterizing 
the strength of pinning, and the function go(t) which is 
defined by the type of pinning. Other physical quantities 
[like T c , H C 2(T), e, etc.] either determine only scaling 
factors or do not appear explicitly in the appropriate 
equations at all. In particular, the anisotropy e is ab- 
sorbed by the definitions of Gi and D and does not enter 
in Eq. Pty . 



B. Analysis of the equation 

Figures [21 and |3| present the melting lines H m (t) ob- 
tained numerically from Eq. (|39|) for different pinning 
strength D. Although the numerical factors A and 
B have remained unknown in the above derivation of 
Eq. I|39l) . the figures show that it is possible to choose 
Ax = A/(2tt) 3 / 4 and B x = B/i/(27r) 3 / 2 so that the cal- 
culated melting lines have the properties observed in ex- 
periments (see Introduction). In particular, we find that 
the upper critical point noticeably shifts under the influ- 
ence of pinning and tends to T c at reasonable values of 
D when A\ and B\ are not too small (A±, \[~B\ ~ 1). 
With increasing D the melting line shifts downwards if 
Bi < A\ , and this shift becomes small when Bi does not 
differ too much from A 2 . Of course, the choice of A\ and 
B\ on the basis of these requirements is not unique. Let 
us now analyze Eq. (|39J) in some limiting cases. 

According to Sec. II, in the temperature region (fl*2"|l 
one has 1 — h m <C 1 for the melting line of the ideal 
vortex lattice. If the parameter D is not too large, this 
property of the melting line remains true for the vortex 
lattice with pinning. But under the condition l — h m <C 1 
the lowest Landau level approximation is valid, and we 
can use the results of Appendix A where the effect of 
quenched disorder on the melting line has been analyzed 
in the framework of an approach similar to that of Sec. II 
C. On the other hand, putting h m w 1 in Eq. (|39|l . we 
arrive at 



-3/2 



Bfi[G(t)} 2 
[I + O- 3 / 2 ] 2 
AG{t)Q~ z / 4 
[I + Q- 3 / 2 ] 1 / 2 



27rc|. 



(43) 
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where Q 3 ! 2 
dHJ], and 

G(t) 



(1 — h m ) 3 / 2 1 Fx(t) [compare with formula 



[Dg (t)] 3 / 2 [Dg (t)] 3 / 2 (l-t 2 )^ 



[FT(t)} 1/2 



(44) 



Equation 143|l shows that the quantity Q on the melting 
line depends on t only via the function G(t). This conclu- 
sion completely agrees with the exact result of Appendix 
A. 

When D = [i.e., G(t) = 0], Eq. g3J) gives the result 
for the ideal vortex lattice, Q~ 3 ^ 2 = 2-kc\, see Eq. ijTT?)) . 
For nonzero D the function G(t) decreases with the tem- 
perature. In the case of 5T C pinning one has G(t) cx t~ x / 2 , 
while G(t) is proportional to (1 — t 2 )/^/t for 51 pinning. 
Hence in the equation for the melting line the role of 
terms associated with pinning increases with decreasing 
t. It follows from Eqs. (|4*3l . (|4*21) that at the upper crit- 
ical point the function G(t) reaches a certain value, and 
condition (|42[1 can be written in the form: 



(27TC 2 



• Bfi 
A 2 



1/2 



(45) 



Interestingly, up to a numerical factor (and the factor [1 — 
h sv (t)} 3 / 2 which we do not take into account here^) this 
condition reduces to Eq. (32) of Ref.0for the intersection 
point of the melting line with the order-disorder line. 
In particular, it follows from Eqs. (|44|l . H45|) that t up , 
similarly to the temperature of the intersection point, 
depends on Gi and D only via the combination D 3 /Gi 1 / 2 
of these parameters. 

The so-called depinning line 9 along which ut ~ £, in- 
tersects the melting line of the ideal lattice, H^(t), at 
1 — t 2 ~ (f 2 /2ttc 2 j ) Gi, i.e., outside the temperature 
region defined by Eq. ltT2)l . In the region 1 — t 2 <C 
(f 2 /2nc 2 j ) Gi, the depinning is essential. In this depin- 
ning region, the terms TAS and E pin in the energy bal- 
ance (|27|) become relatively small, and the melting line 
is close to that of the ideal lattice. 



C. Lindemann criterion 

We now rewrite Eq. i|39[) in the form of the Lindemann 
criterion. As it follows from formulas © and 10, the 
factor F T {t)h]l 2 /(l - h m ) 3/2 in the first term of Eq. ^ 
is simply 2-nu\.ja 2 where ut is described^ 3 - by Eq. JBJ). 
To rewrite the second term in Eq. I|39fl . let us take into 
account the formula for the averaged relative shift of two 
line elements in the vortex lattice separated by a distance 
R < R c transverse to the magnetic field^i 



u 2 {R, 0) 



eafyi + ln(R 2 /a 2 ) + eR/X ab 



L c r 2 



(1 - ft)3/2 



(46) 



[If the small bundle pinning regime occurs, R c < X a b, 
the last term in this formula, eR/X a b, is small and may 
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FIG. 2: The melting line H m (t) for the vortex lattice with 
quenched disorder, from Eq. l Eflfl . for D/c L = 0, 0.4, 0.6, 0.8, 
1, 1.2, 1.4, 1.6 in the case of 51 pinning at Gi = 0.01 (top) 
and at Gi = 0.001 (bottom). Here Ax = A/(2tt) 3/4 = 1. 
The upper line D = is the melting line of the ideal lattice, 
H^(t). The dotted line shows H c2 {t)/H c2 {0) 



be omitted.] At R — R c , one has u(R c ,0) 



r p . Insert- 



ing r 2 = ? 



and the expression for ut in the 



relationship thus obtained, and taking into account the 
definition i|32|l . we arrive at 



R 2 

2 



1 In ( — ] +3^ 



X a b 
1 



h(l - h) 
^D 2 g 2 {t) 



3/2 



(47) 



1/2 



fit 



Gi 



l-t 2 J /l 1 /2(l_/ l )3/2 



Simple manipulations using this formula show that the 
second term in Eq. (|39fl , which tells the role of the pinning 
energy in the energy balance, is 

-1/2 



2-kA x 



1 + ln 



Rf 



eR c 

X a b 



(48) 



where A\ = A/(2t:) 3 ^ 4 . Since the factor containing R c 
can be represented as u(a, 0)/r p [using Eq. (|4*r)|) at R = a 
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FIG. 3: As Fig. H but for 5T C pinning. 



and then at R = R c ], we find one more form for the 
second term in Eq. (jSHJ: 



2nA 1 



r p u(a, 0) 



(49) 



As to the quantity P(t,h) given by Eq. (|4U[) . it is 
expressed as B 1 u 2 (a, 0)/urj, where Bi = Bfi/(2w) 3 / 2 . 
Thus, equation (|3*9"jl is equivalent to the following cri- 
terion: 



iiy — Biu 2 (a, 0) + Air p u(a, 0) = c|a 2 



(50) 



which is true for the melting line in the bundle pinning 
region. Equation (|50|l can be also rewritten in the form: 



u 2 T - B lU 2 (a, 0) + A lU 2 (a, 0) 



\'4 

a 



eR c 

A a 6 



1/2 



c\a 2 . (51) 



Note that although criteria (J2J and J3J) are qualitatively 
close to Eq. (|51|) . they underestimate pinning at large R c . 

Equation [or l|5T|) ] is valid when u\ > Biu 2 (a,0) 
since condition l|42|l for the upper critical point now has 



B lU 2 {a,0). 



(52) 



Taking into account that B± ~ 1, this condition means 
that the upper critical point, as a rule, lies in the bundle 
pinning region. Indeed, if the melting line enters the sin- 
gle vortex pinning region before the upper critical point 



occurs, one finds u 2 -. > B\r 2 = £?i(£ 



at the bound- 



ary of this region. In other words, the melting line can 
intersect the boundary in the depinning region where u\ 
is large. But in this region the melting line practically co- 
incides with that of the ideal lattice, Sec. Ill B; the latter 
can cross the boundary only in the case of 8T C pinning 
at sufficiently large 



IV. PHASE DIAGRAMS 

A. The first and the second scenarios 

To gain some insight into the character of the mergence 
of the melting line with the order-disorder line, let us 
consider more closely the dependence of the free energy 
difference between the liquid and the Bragg glass, AF, 
on the size of the dislocation cell, l±. We shall consider 
only such l± that lie near the minimum possible size in 
the lattice Z™ n (/" nn ~ a) . In evaluating the entropy gain 
TAS, the elastic energy E c \ and the pinning energy E p - m , 
it is necessary to take into account that the modulus C44 
depends on the wave vector k. According to Ref. |25|, one 
has C44(fc ~ l/'±) ~ 044(1/0) • (/j_/a) 2 at l± < \ a b/e, and 
hence oc l\ : see formula (|17|l . Then we can write AF 



in the form (we still neglect logarithmic factors): 
AF cx E el {a)l 2 ± - E pin {a)? ± + T[P{a)l\ 
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(53) 



P(a) denotes P(t,h m ) of Eq. [j39"j> . 

pin (a) are the elastic and the pinning 



where l ± = l±/l^ ln 
and -E e i (a) and E ] 
energies per dislocation cell at l± = l™ m , which have been 
inserted into formula (|27[1 to derive Eq. (|39[) . Equation 
(|53fl is valid when P(a)l 4 L — 1 < 0, otherwise the last term 
in this equation has to be omitted, and one has, 



AF oc E e i(a)Tj_ — E p - m (a)Tj_, 



(54) 



for l\ > 1/P(a). 

As was mentioned above, in the equation (|39|l for 
H m (t) [or equivalently T m {H)] the relative role of terms 
associated with pinning increases with decreasing T m for 
both types of pinning. When F p j n (a) and P(a) are suffi- 
ciently small (i.e., for temperatures T m near T c ), the func- 
tion AF(^j_) is minimum at the lowest possible value l± = 
1 and then increases with increasing l±, Fig. This 
means that melting occurs at l± = 1 as in the ideal lat- 
tice. This conclusion has been already used implicitly in 
deriving Eq. When the temperature T m decreases, 

and hence P(a) and the ratio F p ; n (a) / ' E c \(a) increase, the 
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point Zj_ = [l/P^)] 1 / 4 shifts to smaller l±, and the slope 
of AF in Eq. I|54|) on l\ decreases. Eventually, at the crit- 
ical temperature t up , one arrives at the situation shown 
either in Fig. 0t or in Fig. 0Ji since the energy balance 
for the melting E e \{a) — E p i n (a) = T m [l — P(a)] leads to 
E e i(a) — E p i n (a) = at this point. In the case of Fig. ^Jl 
when the dependence of AF on P\ has a negative curva- 
ture at lj_ — 1, we conclude that this case is necessarily 
preceded by the situation shown in Fig.^Ja, which has to 
occur at some temperature U > t up where P(a) < 1 and 
E p in(a) < -Eei (a). At this U the size of the dislocation 
cell, l± , sharply increases because the absolute minimum 
of AF now occurs at a larger l± > 1. In other words, at 
this temperature we find an intersection of the melting 
line with the order-disorder line. The abrupt increase of 
l± also means that the melting and the order-disorder 
transition cannot form a unified phase transition line. 
Thus, we arrive at the first scenarios^ shown in Fig. ^ 
In the case of Fig. two possibilities exist: Either at 
some temperature U > t up the situation shown in Fig. 2t> 
occurs, and we again arrive at the first scenario, or the 
function AF(l±) — AF(1) remains positive for all l± in 
the temperature interval t up < t < 1. In the latter case 
we find that the mean size of the dislocation cell, l±, can- 
not gradually change in the process of the reduction of 
the entropy gain, and hence one has l± — 1 down to t up . 
In this case, the melting line continuously transforms into 
the order-disorder line, and we arrive at the second sce- 
nario. Moreover, since t up is not a specific temperature 
for Spin (a) or E c \(a), it is quite probable that the result 
l± = 1 remains true also at t < t up (at least in the bundle 
pinning region). Then, at t < t up the order-disorder line 
(which is now a part of the unified phase transition line) 
is determined by the condition, 



E cl (a) = E pin (a). 



(55) 



As it follows from the results of Sec. Ill C, in the bundle 
pinning region condition (|55|) is equivalent to the Lindc- 
mann criterion: 



A lU 2 (a,0) 



1 + ln 



P, 



eR c 



1/2 



2 2 

c L a . 



(56) 



If the order-disorder line enters the single vortex pinning 
region, where R c = a, condition l|55ll yields: 



AiM 2 (a,0) = c\a? 



(57) 



(Note that the ratio ea/X a b is negligible.) In Appendix 
B equations (|56|l and (|57J) are presented in explicit form. 

To calculate the order-disorder line H<n s (t) in the case 
of the first scenario, we used the criterion^ 



u 2 (o,0) 



-LD L 



(58) 



where cld is some new Lindemann constant describing 
the order-disorder transition. It should be emphasized 
that although Eq. l(5"%)) is similar to Eq. (|S*7|) . the Lin- 
demann constant cld in Eq. (|58(l is independent of the 




FIG. 4: Schematic sketch of the function AF(l±), Eqs. t55l . 

for t > U (a), t = ti (b), and t = t up (c or d). The 
parts of the curves where Eqs. 15311 . 1541 fail are indicated as 
dashed lines. The breaks of the curves a and b occur at the 
points l\ — 1/P(a) where Eq. 1531 transforms into Eq. 1541 . 



constant cl defining the melting, while in the case of 
the second scenario Eq. I|57(l leads to a relationship be- 
tween the Lindemann constants defining the melting and 
the order-disorder lines (the ratio of these constants is 
V^4i). Criterion (|58|l both in the single vortex and in the 
bundle pinning^ regions were analyzed in Ref. 0, allow- 
ing for a smoothing of the pinning potential by thermal 
fluctuations. The appropriate equations in the explicit 
form are presented in Appendix B. For definiteness, in 
the analysis below we shall choose the constant c^d for 
the order-disorder line as equal to the constant cl for 
the melting line. 



B. Numerical results 

Although in this paper we have not analyzed the order- 
disorder line in detail, criteria l(5T)|) . (|5Tjl . and (|55|l enable 
one to evaluate the location of this line in the T - H 
plane for the first and the second scenarios. In Figs.[S]- 
|S]we compare the T-H phase diagrams of type-II super- 
conductors with different types of pinning for these two 
scenarios. In the construction of the figures, we use the 
values of the Ginzburg number Gi ~ 0.01, 0.0001 which 
are typical for high-T c superconductors, and we take into 
account the factor containing [1 — h sv (t)] that has been 
omitted in Sec. Hl4£ The complete set of the appropriate 
equations is given in Appendix B. 

In Figs. [S] - |H1 we choose the constants Ai and B\ so 
that the following two requirements are satisfied: First, 
the properties of H m (t) observed in experiments are re- 
produced, namely, with increasing D the upper critical 
point clearly shifts, while the downward shift of the melt- 
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ing line is small, see Sec. Ill B. Second, the upper crit- 
ical point t up of the melting line coincides with the in- 
tersection point ti of this line with the order-disorder 
line (calculated within the first scenario) at D — 0.7 for 
the case of 51 pinning. The latter requirement is due to 
the following considerations: According to Ref. 0, 0, 
the coincidence of U and t np is observed in overdoped 
YBa 2 Cu30 y crystals (y > 6.92) for which the upper crit- 
ical point lies at sufficiently large magnetic fields, and so 
the case h m ~ 1 appears to occur there. As it was men- 
tioned in Sec. Ill B, in this case U and t up depend on Gi 
and D only via a single combination of these parameters, 
D 3 /Gi 1 / 2 . Therefore, one may expect that if the coin- 
cidence of the upper critical point with the intersection 
point occurs, it practically will not depend on the spe- 
cific choice of D or Gi. In other words, the coincidence 
will approximately occur for different D, Gi, and types of 
pinning as long as fijRil, Although hi is not too close to 
unity in Figs. EJ-EI (top), the data of these figures support 
this statement. As one might expect, the region of the 
coincidence is especially wide in D for small Gi, Fig. [5] 
Thus, we have introduced the second requirement here in 
order to fit the phase diagrams calculated in the frame- 
work of the first scenario to the experimental situation 
observed in the overdoped YBa2Cu30 y crystals. 

The presented figures show that when the intersection 
point is close to the upper critical point, both scenar- 
ios lead to qualitatively similar phase diagrams. How- 
ever, when the strength of pinning D increases, and the 
intersection point shifts toward T c , the coincidence of 
the points fails especially for ST C pinning. In this case 
the first scenario leads to a noticeable extension of the 
melting line beyond the intersection point. Thus, if the 
first scenario really occurs, this result possibly explains 
the experimental findingsii*i2iiii^ii^ for optimally doped 
YBa2Cu30j, crystals (y « 6.92) in which the upper criti- 
cal point lies at larger magnetic field than the intersection 
point. Note also that in the case of ST C pinning, the melt- 
ing line enters the single vortex pinning region near T c at 
v = (27r) 3 / 2 D z /Gi 1 ! 2 of the order of several units4*i^ 

In Figs. - 0D we also show the boundary of the sin- 
gle vortex pinning region, H sv {t). It should be noted 
that apart from the well-known lower single vortex pin- 
ning region^ an upper region exists where this pinning 
occurs. The upper region of single-vortex pinning was 
discussed by Larkin and Ovchinnikov 47 in the context of 
the origin of the peak effect in low-T c superconductors. 
Without account of thermal fluctuations, this region ad- 
joins to the H C 2(t) line and is caused by the softening 
of the vortex lattice near H C 2(t). However, this softening 
also leads to an increase of Ut, which reduces the strength 
of pinning. As a result of these two opposite tendencies, 
the upper region does not extend to T c and has the shape 
of a "tongue" . 4 Interestingly, when the strength of pin- 
ning increases, the upper and lower single- vortex pinning 
regions can merge at low temperatures. This merging 
occurs for D > (47r) -1 / 2 ; this value is independent of 
Gi and the type of pinning since the merging starts at 
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FIG. 5: The phase diagram in the case of SI pinning, go(t) = 
{l-t 2 ) 1/2 . Here Ai = 0.66, B = O.&A?, c L = 0.25, Gi = 0.01, 
and D/cl = 0.8 (top) or D/cl = 1.3 (bottom). The melting 
line H m (t) and the order-disorder line H^ a (t) for the first 
scenario are shown by solid lines, while the thick dotted line 
gives the order-disorder line H^ a (t) for the second scenario. 
The dashed line depicts the boundary of the single vortex 
pinning region, H av (t), Eq. (B2) [the thin-dashed line shows 
Eq. (Bl)]. The dotted line is H c2 (t)/H c2 (Q) = 1 - t 2 . The 
upper critical point (T up , -Hup) is marked by a dot and the 
intersection point (Ti, Hi) by a circle. 



T = 0. 

The authors of Ref. argued that only the second sce- 
nario can explain the decrease of Hdi S with the temper- 
ature that was observed in their experiments. However, 
the presented figures show that the order-disorder line 
found from criterion (|58|l (the first scenario) can decrease 
or increase with temperature and even can be nonmono- 
tonic. This depends on whether the line is in the single 
vortex pinning region or in the bundle pinning region 
and also on the Ginzburg number. It should be noted 
that in contrast to Eq. © we take into account only the 
displacement caused by the quenched disorder and do 
not include ut explicitly in the equation for the order- 
disorder line, compare Eq. (|58|l with Eq. ©. However, 
the thermal depinning (which depends on ut), the soft- 
ening of the elastic moduli, and the possibility that the 
order-disorder line lies not only in the single vortex pin- 
ning region but also outside it, already produce the de- 
picted variety of shapes of Hdi S {t). Hence, even in the 
framework of the first scenario the presented results can 
explain the facl£ that the order-disorder lines observed 
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FIG. 6: As Fig. but for Gi = 0.0001. 



FIG. 8: As Fig. □ but for D/c L = 1.15. The lower panel 
shows the same phase diagram near T c at enlarged scale. 
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FIG. 7: As Fig. but for the case of 5T C pinning, g {t) 
(1 - t 2 y 1/& . Top: D/c L = 0.8; bottom: D/c L = 1.12. 



in experiments have various shapes. 

Finally, we briefly discuss the case of low-T c super- 
conductors, which have a very small Ginzburg number. 
In this case one has D 3 /Gi 1 ^ 2 1 even for weak pin- 
ning strength D. Since the temperatures Ti and T up are 
mainly determined by this ratio of the parameters and 
increase when the ratio increases, we find that in low- 
Tc superconductors these temperatures practically coin- 
cide with T c , and only the order-disorder line can be 
observed^ Interestingly, in this case the order-disorder 
lines for the first and for the second scenarios have the 
same functional dependences on t and on the parameters 
D, Gi [this follows from equations l|B5|) and l|B8|) ]. and 
moreover, they practically coincide with each other when 
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APPENDIX A: EFFECT OF DISORDER ON THE 
MELTING LINE: ANALYSIS BEYOND 
PERTURBATION THEORY. 

Assuming 1 — t -C 1, we shall consider the partition 
function of the vortex system as a functional integral with 
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the Ginzburg-Landau Hamiltoniani2S^4 In dimensionless 
units this Hamiltonian has the form: 



H GL 



dr 



|9^| 2 + l(-»V + A)^| 2 + r|^| 2 + 

2 i ' 2 nGi 1 ' 2 m' l \ , (Al) 



where ip is the order parameter; r = t — 1; the coor- 
dinates x and y are measured in units of £o> an d z m 
units of e£o! Co is the zero-temperature coherence length 
in the Ginzburg-Landau theory; the magnetic field is b = 
2it£ i qH/$> an d A is its vector potential; V = (d x ,d y ). 
Consider the melting line h m (t) in the temperature region 
described by Eq. l(T2)l in which 1 — h m <C 1, and hence 
the lowest Landau level approximation is valid (see Sec. II 
C). In this approximation the second term of Eq. (|A1|) 
reduces to 



(A2) 



Pinning is introduced into the dimensional Ginzburg - 
Landau Hamiltonian either via spatial disorder in the 
transition temperature T c + 6T c (r) (5T C pinning) or by 
spatial variation of the effective mass m + 5m(r) describ- 
ing disorder in the mean free path I of quasiparticles (51 
pinning) f& Here m is the effective mass in the x-y plane 
(since the magnetic field is along the z axis). Thus, with 
quenched disorder in the vortex lattice, one should add 
to the dimensionless Hamiltonian (|A1|) the term 



where 



<p(r) = 



6T c (r) 



in the case of 5T C pinning, and 



^ (r) = M^) 6Rj Mr) |r| 



(A3) 



(A4) 



(A5) 



in the case of 61 pinning. In Eq. I|A5() we have used the 
relation 



Mh 



(A6) 



that follows from the definitions of b and h and put h = 1 
since 1 — h -C 1 in the lowest Landau level approxima- 
tion. For pinning by point defects, it is assumed in the 
collective pinning theory 9 that disorder in ST c (r) and in 
<5m (r) is short scale and described by a Gaussian distri- 
bution with zero mean value, (<5T c (r)) = (Sm(r)) = 0, 
and with the correlation function^ 



(6T c (r)6T c (r')) 



rp2 
± C 

for ST C pinning, and 



27rL> 3 (5(r - r') 



(5m(r)5m(r'))_ = 307r g3 _ . 



(A7) 



(A8) 



for 81 pinning. Here (. . .) means disorder averaging. Note 
that in agreement with Ref. and with Sec. Ill A of this 
paper, equations IjAlfl - i|A8(l again show that the phase 
diagrams of type-II superconductors with point defects 
depend only on D and Gi. 

To proceed further, let us rescale the coordinates and 
the order parameter ip similarly to Ref. l35t 



x = Vbx, 
V = Vby, 



1 = zGi^Htb) 1 '*, 
4> 2 = ^Gz^t^b- 2 ^, 



(A9) 
(A10) 
(All) 
(A12) 



where x, y, z, are the new coordinates and order pa- 
rameter. Beside this, to agree with the notation used in 
the main text of this paper, we put r = t 2 ~ 1 below. 
Then, the Hamiltonian (|A1|) transforms into: 



H, 



GL 



T 



dr 



I^| 2 -QM 2 + 2 3 / 2 tt|M (A13) 



where Q is defined by Eq. (|21|l . In the case of the ideal 
lattice, Eq. (|23|) for h m follows from this Hamiltonian, 
see Sec. II C. In a lattice with quenched disorder, the 
additional term from (|A3|) in the Hamiltonian has the 
form: 



mm 2 , 



with (</?(?)) = and 



27rD 3 

in the case of ST C pinning, and with 

,- rw -m 30^P 3 (l-t 2 ) 2 „ 

WW )> = 7Gi i/2 t S(t-t) 



(A14) 



(A15) 



(A16) 



in the case of 51 pinning. 

Up to numerical factors, the right hand sides of 
Eqs. (|A15j) and (|A16jl coincide with the function [G(t)] 2 
defined by Eq. I|44fl . Then, the free energies of the vortex 
liquid, i^i q , and the vortex lattice, Fi at , are determined 
by Q and G(t), and along the melting line h m (t), one 
has: 



■Fiiq(Q, G(t)) = FutiQ, G(t)), 



(A17) 



Thus, on the melting line, the quantity Q is some function 
of G(t). Note that the functional form of Eq. Ij43(l agrees 
with this conclusion, see Sec. Ill B. 



APPENDIX B: EQUATIONS FOR 
CALCULATION OF THE PHASE DIAGRAMS 

In section III where we considered the melting line in 
the bundle pinning region, the normalization factor con- 
taining [1 — h sv (t)\ was omitted^ The origin of this fac- 
tor is due to the additional power of 1 — h in the shear 
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modulus C66 as compared with the tilt modulus C44, see 
Eqs. O. In the single vortex pinning region where the 
shear modulus does not play any role, the additional 
power should not manifest itself in physical properties. 
Therefore, if a physical quantity in the bundle pinning 
region contains a factor (l — h) n caused by this additional 
power, it is necessary to introduce a normalization factor 
[1 — h sv (t)]~ n in this quantity to provide its continuity at 
the boundary of the single vortex pinning region h sv (t) . 
(Here n is some power.) Although the effect of the fac- 
tor [1 — h sv (t)]~ n on the melting line H m {t) and on the 
order-disorder line H dis (t) in the bundle pinning regime 
is small, it essentially influences the boundary of the sin- 
gle vortex pinning region, H sv (t), and the line H d i S (t) in- 
side this region. In this Appendix B, taking into account 
this factor^ we compile the complete set of equations 
used for the construction of Figs. 0-03 

The boundary of the single vortex pinning region 
h sv (t) = H sv (t)/H c2 (t) is described by equation (19) of 
Ref. 4: 



h\' v \t) 



F T (t) 



[l-h sv (t)}*/< 



= {2^) 1 ' 2 Dg {t) , (Bl) 



where we have used the notation i|4ip. The boundary 
of the upper single vortex pinning region, discussed in 
Sec. IV B, can be obtained from the equation: 4 



1 



1 - h sv (t) 



F T (t) 
(1 - hf' 2 



= 2tt (£><?[,(*)) , (B2) 



which in addition reproduces the root of Eq. (jBljl . i.e., 
it yields the entire boundary of the single vortex pinning 
regions. Inside this upper region the vortex lattice is in 
a state where R c = a and u(a, 0) = r p , i.e., a border- 
line state between the single vortex pinning and bundle 
pinning regimes occurs there. 

Equation (|39l) for the melting line in the bundle pin- 
ning region is rewritten as follows: 



F T (t)h 



1/2 



[1-P(t,h m )] + 



(1 - M 3/2 

A[Dg {t)f' 2 hU 4 [l~h sv (t)] 3/i 
(l-h m )V 2 +F T (t)h- 1/211 1 



2nc 2 L . 



(B3) 



with 



P(t, h m ) = 



B[Dg (t)] 3 (l-t 2 ) l / 2 [l-h sv (t)] 



tGi 1 / 2 



1/2 



F T {t) 



1 2 



(1 - M 3/2 



(B4) 



Equation (|B3() is valid in the interval t up < t < 1 where 
the temperature t up defines the position of the upper crit- 
ical point of the melting line and is given by the condition 



(I42[l . Note that in this paper we consider situations when 
the upper critical point lies in the bundle pinning region. 

In the framework of the second scenario the equation 
for the order-disorder transition line (which is the con- 
tinuation of the melting line to t < t up ) in the bundle 
pinning region follows from criterion (|56[) as: 



A[Dg (t)}^ 2 hU 4 [l-h s 



(1 - M 3/2 + F T {t) h] 



^ = 2^1, (B5) 



-1/2 



i.e., the first term in Eq. i|B3|l disappears. If this order- 
disorder line enters the single vortex pinning regions, it 
is described by equations that result from criterion 15711 . 
This criterion reads in explicit form 



A 1 h 3 J v 5 (t)h 2 i 5 



1 



F T {t) 



hli 2 (l-h m )3/ 2 



= 2-kcI (B6) 



for the order-disorder line in the lower single vortex pin- 
ning region, while in the upper single vortex pinning re- 
gion one finds: 



MF T (t) 



,1/2 



(1 - M 3 / 2 



Ai h m = 2nc 2 L . 



(B7) 



Here we have used that u(a, 0) = r p in the upper region. 

In the framework of the first scenario, equations for 
the order-disorder line follow from condition (|58|l . 4 With 
the use of our approximation /1 = 2, equation (24) of 
Ref. for the order-disorder line Hdi S (t) in the bundle 
pinning region can be written in the form: 



i 1/2 

l dis 



(1 



s.3/2 



F T (t)K±{t), 



where h dis = H dls (t) / H c2 (t) , 

K ± {t) = G 2 -l±[(G 2 -l) 2 -l} 1/2 



and 



G x = 



^l 4 G{t)\\-h sv (t)fl 4 



2V4 



CL 



(B8) 



(B9) 



(BIO) 



with G(t) from Eq. glj. The order-disorder line H dis {t) 
in the single vortex pinning region is given by Eq. (22) 
ofRef.S 



F T (t) 



{i-h dls f/ 2 [h sv {t)Yi 2 



h sv {t) 



2ttc1 



5/2 



3/2 



(Bll) 



Note that apart from the factor A\ Eq. (|B11|) is equiva- 
lent to Eq. (|B6|) . 
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